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1. Introduction 

Early in the 21st century space infrared observatories such as ASTRO-F,-'^ Space Infrared Telescope Facility 
(SIRTF),^ Far Infrared and Submillimeter Telescope (FIRST), and Space Infrared Telescope for Cosmology 
and Astrophysics (SPICA)'* are planned to study the formation of planets, stars, and galaxies. The improve- 
ment in sensitivity is so great that the source detection limits of these space observatories are expected to 
be set by source confusion, especially in the far infrared. To further improve the source detection limits, 
individual sources need to be resolved by an infrared interferometer. Therefore it is important to quantify 
the sensitivity of a space infrared interferometer for the purpose of synthesis imaging. 

Following the observatories, space infrared interferometers such as the Terrestrial Planet Finder (TPF)^ 
and Darwin^ are considered, with an emphasis on the detection of terrestrial planets around nearby stars. 
TPF and Darwin will consist of several radiation cooled apertures and will be capable of synthesis imaging 
as well as planet detection by nulling interfcrometry. The detection limits of these interferometers for general 
synthesis imaging are of great interest. Among many applications of synthesis imaging with a space infrared 
interferometer, we are particularly interested in deep galaxy count by which the history of galaxy formation 
can be studied. We discuss how a TPF-like interferometer operated in the far infrared can improve the 
number count data of high redshift galaxies. 

The beam combination geometry is a major issue in studying the sensitivity of an interferometer. One 
extreme is the "C2 interferometer, in which n beams are divided into n(n— 1) subbeams and they arc pairwise 
combined at "C2 — n{n—l)/2 detectors. There is one detector for each baseline. The other extreme is the "C„ 
interferometer in which all the n beams are combined at "C„ = 1 detector. Prasad and Kulkarni^ (henceforth 
PK89) have rigorously analyzed the SNRs of the "C2 and "C„ interferometers for spectrally narrow light in 
the presence of source shot noise and concluded that the sensitivities of both interferometers are essentially 
the same apart from a numerical factor. The results by PK89 are applicable to optical interfcrometry, but 
not to infrared interfcrometry because of the presence of thermal background and detector read noise. The 
beam combination geometry can have an impact on the sensitivity in the infrared for instance because of 
the number of detector pixels which affects the total read noise. The beam combination geometry may also 
have a major impact on the architecture of the space interferometer. 
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Baseline redundancy is a factor that affects the SNR of a synthesis image for a given number of apertures. 
Since the number of apertures deployed or flown separately in space is likely to be limited, it is important to 
gain uv coverage by using nonredundant baselines. Throughout this paper, we only consider interferometers 
with nonredundant baselines. 

There are two types of noise which appear in aperture synthesis, one of which is removable while the 
other is intrinsic and not removable. Sampling noise is potentially introduced by the incompleteness of uv 
coverage. However, deconvolution techniques such as CLEAN and Maximum Entropy Method (MEM),^ when 
applied to radio interferometric data, appear to compensate for sampling noise. We thus regard sampling 
noise as removable noise. On the other hand, noise resulting from the photodetection process is not removable. 
The SNR of a dirty image to which deconvolution has not been applied is limited by photodetection noise. 
Here we derive the SNR of a dirty image as determined by source shot noise, thermal background shot noise, 
and detector read noise, and equate it to the SNR of the synthesized image. 

There are two different methods by which a synthesized image can be constructed from the visibility 
data: inversion without total counts and true inversion. In the former, the zero frequency phasor is neglected, 
implying that the total number of photons in the synthesized image is zero. Despite this seemingly unattrac- 
tive feature, this is the standard method in radio astronomy. In the latter method, the van Cittert-Zernike 
theorem is strictly applied and the zero frequency phasor as well as the n{n — l)/2 complex phasors is 
used. The image produced by this technique has the desirable property of nonnegativity. In the presence of 
additive thermal background which overwhelms the signal, the magnitude of the zero frequency phasor is 
much greater than those of other phasors. The zero frequency phasor determines the DC offset level due to 
thermal background in the synthesized image. Not only is this background of little astronomical interest, 
it is also a major source of Poisson fluctuations. Therefore in this paper we prefer the former method and 
consider only the inversion without total counts as in radio astronomy. 

We analyze the SNRs of space interferometers in three steps. First, we formulate the SNRs of "C2 
and "C„ interferometers in the presence of source shot noise, thermal background shot noise and detector 
read noise in §2 and §3 respectively. These formulae cover optical and infrared interferometers. Second, we 
consider spectrally broad light for which fringe dispersion becomes an issue. Fringes of the "C2 interferom- 
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cter can always be dispersed while those of the "C„ interferometer can be dispersed only in the case of a 
one dimensional baseline configuration. We discuss the tradeoff between the "C2 and "C„ interferometers 
including the dimensionality of the baseline configuration in §4. Third, we apply the results of our analysis 
to planned interferometers. We present the sensitivities of a TPF-like interferometer and a Darwin- like one 
in §5. Finally in §6 we discuss the capability of the TPF-like interferometer from the point of view of deep 
galaxy count. 

2. "C2 Interferometer 

A. Fringe Phasor Estimator 

Let there be n identical principal apertures from which we derive n main beams. Each beam is divided into 
n — 1 subbeams by the use of beam splitters. The resulting n{n — 1) subbeams are combined pairwise on 
rib = n{n ~ l)/2 detectors where rtf, is also the number of baselines. 

The fringe pattern is formed at the focal plane of each detector. For spectrally narrow light, a one- 
dimensional detector is sufficient (Fig. 1). For spectrally broad light, the fringes can be dispersed in the cross 
fringe direction, in which case a two-dimensional detector is used. For a focal plane interferometer of this 
type, the influence of source shot noise on fringe phasor estimation has been fully formulated^ Let the 
interferometer be illuminated by a source upon a spatially smooth background. The intensity pattern on the 
r'^ detector (r"^ baseline) is given by 

< /r(x) >= 2 < /o > 1 + 7r C0s(kX • Br/d + 0r) , (1) 

where < /o > is the average intensity in each subbeam at the detector, is the baseline vector, 
K ~ 27r/A is the wave number, and d is the distance between the aperture plane and the detector, x is the 
spatial vector in the detector plane, and 7rexp(i</)r) is the complex visibility function at the baseline vector 
B^. In deriving (1), we have assumed that the incident light is spectrally narrow so that the fringe visibility 
depends only on the spatial correlation of the field. 

Because of the presence of background illumination. 
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lo^I^+l'o, (2) 
where Iq and /q are source and background intensities respectively. Let 7^ be the fringe visibihty of the 
source in the absence of the background, then 



> = -Wh^^r- (3) 

The complex visibility function 7,^ exp{i(f>r) is determined by the uv coordinates of the r^^ baseline {u, v) 
and the source brightness distribution on the sky S{x,y) by a Fourier transform relation: 



Pxudrh ) - / Sjx, y) eM-^M^x + vy)}dxdy 
7,expiz0,j- JS{x,y)dxdy ' ^' 

In an effort to reduce the clutter in the equations we henceforth drop the vector notation, but bear in 

mind that spatial frequencies, pixel locations, etc. are really vectors. The photoelectron detection theory^ 

takes into account the discrete nature of both photons and detector pixels. The average photoelectron count 

< kr{p) > at the pixel location specified by the integer index p of the detector is proportional to < I{x) >: 

< kr{p) >~ 2 < Ko > 1 + 7r COs{pLJr + (pr) ■ (5) 

Here, < ... > denotes averaging over the photoelectron-detection process. K^) is the discrete version of 
/q. The product pLDr is understood to be the scalar product of the pixel position vector p and the spatial 
frequency ujr expressed in inverse pixel units. 

Let < C > be the average number of photoelectrons detected by the entire interferometer in one 
integration period, and let 2< > be the average number of photoelectrons per detector of each baseline 
per integration time. Clearly then, < C >= 2 < N > nt,, and thus < N >=< C > /{n{n — 1)}. According to 
(5), the average number of photoelectrons per detector is equal to 2 < Kq > P, and thus < Kq > P =< N >, 
where P is the number of detector pixels. 

Due to read noise (p) , each measurement of the pixel p yields Ir (p) given by 

lr{p) = kr{p) +rr{p), (6) 
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where Triji) is a Gaussian random number with zero mean and a standard deviation a. 
The fringe phasor for the basehne r is operationally defined as 



Vr ^'^Ir (P) CXp{-ipUJr ) ■ (7) 

p=l 

Walkup and Goodman have shown that the quantity 

p 

Zr ~ Kip) exp(— ipwr) (8) 

p=l 

is an optimum estimator of the actual fringe phasor under shot noise limited condition in the presence 
of constant background.^ Since 

p p 

< rr{p) e-K.p{—ipujr) >~ ^ < Trip) > cxp(— ipcj^) ~ 0, (9) 

p—1 p—1 

Ur is an optimum estimator of the fringe phasor in the presence of detector read noise. 
The average of the phasor over many coherent integration interval is given by 

l^r =< 2/r > = < > 7rCXp(i0r). (10) 

B. Inversion 

The synthesized image is the real portion of the Fourier transform of the spatial coherence function: 



ii{q) = Re [^y^ yr cxp(+Wrg) 

= '^^[Re{yr) coii{uJrq) — Im{yr)sm{ujrq)]. (11) 

r 

Index q refers to pixels in the synthesized image; in particular q ranges from — Q/2 to +Q/2, and q = Q 
refers to the central pixel in the synthesized image. The variable ii{q) refers to the image obtained from one 
set of visibility data. Note that the sense of Fourier transform used in (11) is consistent with the definition 
of yr- The mean image Ii{q) is the average of ii{q) and is given by 
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rib 

Ii (q) = Re ^ Yr cxp(+iuJrq) 



(12) 



which by virtue of (10), can be simphfied to yield 



Ii{q) = < N > ^ 7r cos{ojrq + (j)r) 

T 

= <N>Y.-^Y,COs{Urq + cl)r). (13) 

The image Ii (q) is referred to as the dirty image in the parlance of radio astronomy. The dirty image 
is the convolution of the true image and the Fourier transform of the spatial frequency sampling function or 
the dirty beam. A synthesized image can be obtained from the dirty image by deconvolution. 

The variance l^[ii((7)] in the synthesized image ii{q) is given by 



V[ti{q)] = < h{qy > ~ < h{q) >^ 
rib rib 

= ^^{[< Re{Vr)Re{ys) > - < Reivr) >< Reivs) >] cos(tJrg) cos(wsg) 
-[< Re{yr)Im{ys) > - < Re{yr) >< Im{ys) >]coa{uJrq) sm{ujsq) 
-[< Im{yr)Re{ys) > - < Im{yr) X Reijjs) >]sv[i{ijjrq) cosluJsq) 
+ [< Im{yr)Im{ys) > — < Im{yr) >< Im{ys) >] sin(wr(?) sin(ti;s(7)|. (14) 

The problem of estimating the variance of the image reduces to that of estimating three types of 
covariance term: cov[Re{yr), Re{ys)], cov[Re{yr), Im{ys)], and cov[Im{yr), Im{ys)]. 



cov[Re{yr),Re{ys)] = [< i?e(y^)i?e(y,) > - < Re{yr) >< Re{ys) >] 

p P 

= [< ^r{p)lsij>') > - < lr{p) >< ls{p') >] COs(w,p) COs(w,.p'). (15) 

p—1 p' — l 

where 



<lr{p)ls{p') > ~ <lr{p) ><ls{p') > = <K{p)ks{p') > - <K{p) ><ks{p') > 
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+ < kr{p)r,{p') > - <kr{p) >< r,{p') > 
+ < rr{p)k,{p') > - < rr{p) >< ks{p') > 
+ < rr{p)rs{p') > - <rr{p) >< r,(p') > . (16) 

There is no correlation of shot noise between different detectors or between different pixels of the same 
detector. The correlation remains only for the same pixel. There is no correlation of shot noise and read 
noise even for the same detector pixel. And there is no correlation of read noise between different detectors 
or between different pixels of the same detector. The correlation exists only for the same pixel. Therefore 

< K{p)ks{p') > - < fcr(p) >< K{P') >= SrsSpp' < kr{p) >, (17) 

< kr{p)rs{p') > - < kr{p) >< rs{p') >=< rr{p)ks{p') > - < rr{p) >< ks{p') 0, (18) 

and 

< rr[p)rs{p') > - < Trip) >< r,{p') >= SrsSpp'a^. (19) 

Therefore 

< lr{p)h{p') > - < lr{p) >< ls{p') >= 5rs5pp'{< kr{p) > +Cr^), (20) 

and 



cov[Re{yr),Re{ys)] = V[Re{yr)]drs 

= 6rs^^{< kr{p) > +a^)cOs{iUrP)^ 



S,,{<N>+—). (21) 



Similarly, one can show that 



cov[Re{yr), Im{ys)] = cov[Im{yr), Re{ys)] = 0, (22) 



8 



and 



cov[Im{yr), Im{ys)] = Srs{< N > +— ). (23) 
Substituting the nonzero covariance elements (21) and (23) into (14), we obtain 



V[iliq)] = ^|(<^>+ — )cos2(c^,,g) + (<7V>+ — )sin2(cj,q)| 



E«^>+^) 



/ ,r <C'> Pf^^ 

= nfc(<7V>+ — ) = ^— +nb — . (24) 

In (24), the term < C > /2 is due to shot noise and the term nf,P(T^/2 is due to read noise. Now we 
consider the specific case of a point source (7,^ = 1) at the phase center {(pr ~ 0) for which /i(0) = 
Since the source is located at the phase center, the SNR of the central pixel in the image is indicative of the 
SNR in the image: 



/i(0) _{<C> /2){^/m+I>^)} ^25) 



3. "C„ Interferometer 

In an "C,i interferometer, all n beams interfere on a single detector and nf, fringes are superposed. Both the 
baseline configuration and the detector can be either one dimensional or two dimensional. As an example, a 
two dimensional ^C^ interferometer is schematically presented in Fig. 2. A special case is the combination of 
a one dimensional baseline configuration and a two dimensional detector for which the superposed fringes are 
dispersed in the cross fringe direction so that spectrally broad light can be used without bandwidth smearing. 
We call an "C„ interferometer of this type an "C'j interferometer. A '^C'^ interferometer is schematically shown 
in Fig. 3. Application of (7) with different spatial frequencies results in extraction of nt fringe phasors and 
the image may be synthesized in the usual fashion. 
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A. Fringe phasor estimator 

Let the interferometer be composed of n identical apertures and let it be illuminated by a source and 
a spatially smooth background. The classical intensity distribution of the interference pattern for the n 
apertures has the average value 

</(x)>=</o> n + 2j2lghCos{KX-Bgh/d + , (26) 

g<h 

where the various symbols have meanings similar to those in (1). gh denotes the baseline gh correspond- 
ing to the apertures g and h. Let < fc(p) > denote the photoelectron coimt distribution due to < /(x) >. 
As in §2, we discontinue the vector notation, assume that the total number of pixels is P, and note that 
< k{p) > is proportional to < I{x) >: 

< >=< Qo > n + ^'-fghCOs{pujgh + (l>gh} ■ (27) 

g<h 

Here < Qo > has approximately the same meaning as < Kq > in §2. However, since there is no beam 
splitting, < Qq (n — 1) < A'^o >. Due to read noise r(p), each measurement of the pixel p yields l(j>) 
counts: 



l{p)^k{p)+r{p), (28) 

where r{p) is a Gaussian random number with zero mean and a standard deviation cr. 

In a typical setup, one can imagine integrating on a two dimensional detector for a period equal to 
the coherent integration interval. A Fourier transform of the resulting image yields 2nh peaks, pairwise 
symmetrical about the origin, which can be identified with the nf, fringe phasors corresponding to the ni, 
spatial-frequency components. Let ygh be the fringe phasor estimator corresponding to the baseline gh, then 

p 

ygi^ = X! ^(P) exp{-ipujgh)- (29) 

Here ygh is an optimum estimator of the fringe phasor for the baseline gh. The mean phasor on the ij 
baseline is given by 
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P g<h 



^ < r{p) > exp(-ipwy) 



(30) 



Since the baselines are nonredundant, oj^ 7^ ^gh, unless (ij) and (gh) refer to the same baseline. 
Therefore only the g ^ i,h ~ j terms remains in (30): 



Yi,j = =< M > jij exp(i0y ), 



where we define < M P < Qo >=< C > /n. 



B. Inversion 



The dirty image, 12(5) is given by 



ij) cos{quJij) — Im{yij) sin(gaji 



l<3 



Using (31), we find the mean synthesized image to be 



h (<?) = Re Yij exp(+i9w,j ) 



= < A/ > ^ 7y- cos(gWy- + 



i<j 



(31) 



(32) 



(33) 



The derivation of the variance V[i2{q)] in the synthesized image i2{<l) is somewhat lengthy and given 
in Appendix A. The result is: 



VMq)] 



n< M > +Pa^ 



-rib- 



< Ad > (n — 2) jij cos{qLLJi 



nb+ < M > " 7^^- cos(gw,j 
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(34) 



For a point source at the phase center, q ~ 0, jlj — 1, 0^ = and the SNR of the synthesized image is 



h (0) _ <M> m/m + ^^^^ 

(35) 



VVWrn v/(" < M > +Pa2)/2+ <M>{n- 2)m/m + I^)} 

<C> + I',)}^{n-l)/n 



^<C> +Pa^+ <C> + im^in ~ 2)/n} 

4. Tradeoff among "C2, "C„ and "C^ Interferometers 

We now compare the SNRs of the "C2, "C„ and "C^ interferometers in the case of a point source at the 
phase center. We are interested in the detection hmits of these three interferometers and we consider the 
situation in which /q <C /q. The total photoelectron count < C > is mostly due to thermal background. In 
the hmit of /g < /g, the SNR of the "C2 interferometer, SNRi, is given by 

SNR, ^ i<C.>m{Wm (3g) 
v/< Ci > /2 + nbPia2/2 

where < Ci > is the total photoelectron count and Pi is the total number of pixels per detector. The 
SNR of the "C„ interferometer, SNR2 is given by 



SNR. = l^l2im^I}I^, (37) 
V< C2 > +P2a' 

where < C2 > is the total photoelectron count and P. is the total number of detector pixels. 
The SNR of the "C; interferometer, SNR'. is given by 



SNR', . < > {/oV/onvF^ (33) 
v/< > +Pia^ 

where < C2 > is the total photoelectron count and P. is the total number of detector pixels on which 
superposed fringes are dispersed. 

< Ci >, < C2 >, and < C. > depend on the spectral bandwidths, the optical throughputs and detector 
quantum efficiencies of the interferometers. Here for simplicity we assume that the optical throughputs 
and detector quantum efficiencies are the same, though in practice the throughputs may differ due to the 
difference in optical configurations. Pi, P2 and P. are determined by the fringe sampling requirements. 
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Expressions (25) and (35) were derived for spectrally narrow light. In practice this corresponds to the 
observations of emission lines. For bright objects, it is always advantageous to disperse fringes and reconstruct 
images at different wavelengths. Especially for emission line objects, spectral dispersion is essential. In this 
case, (36) and (38) hold for the fringe row corresponding to each spectral bin. For faint objects with broad 
spectral energy distributions, spectral images formed from individual dispersed fringes need to be coadded 
to increase detection sensitivities. Below we confirm that (36) and (38) hold also for coadded images of the 
faint objects. 

We first consider the case of the "C2 interferometer. Let be the number of pixels in the fringe 
direction and Py be that in the cross fringe direction. Therefore Pi = PxPy We assume that the photon 
spectrum within the bandwidth is not steep and the number of photons in one detector row is given by 
< ci >=< Ci > /Py. Then the SNR of an image synthesized from a set of fringes in detector rows 
corresponding to 1/ Py of the spectral bandwidth is given by 

snri = i<^^^>mm'^} . (39) 

We further assume the total bandwidth is narrow enough so that the source structure does not change 
within it. In this case, Py images can be coadded to form a final image. The SNR of the final image is given 

by 



SNRi = y/P^jsnn 

(< ci > mm/i^} 



y/< ci > /2 + nbPxCr^/2 

(< Ci > /2){/g//„n 
v/< Ci > /2 + nbP,Pyay2 

(< ci > /2)m/m 



(40) 



V< Ci > /2 + nbPiay2 

Therefore (36) holds for spectrally broad light. Similarly one can show that (38) holds for spectrally 
broad light. 
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A. Background Limit For Spectrally Narrow Light 

For the "C2 interferometer, the condition for background hmited observations is < Ci >^ ni,Pi<j^. In this 
case, the SNR is approximated by 



SNRi 



<Ci> /g 
2 T'' ' 



(41) 



For the "C„ interferometer, the condition for background hmit is < C2 >^ P-zc^- In this case, the SNR 
is approximated by 



bmce 1 < ^2(n- l)/n < V2, for n > 3, the "C„ interferometer is preferred in this case. However they 
differ by at most a factor -\/2 for a large n. There is no difference between the "C„ interferometer and the 
"C^ interferometer for spectraUy narrow hght. 

B. Background Limit for Spectrally Broad Light 

The fractional bandwidth of the "C2 interferometer, Aiz/iy is limited only by the wavelength dependence of 
the astronomical source structure. Each fringe can be spectrally dispersed on a two dimensional detector to 
avoid bandwidth smearing. In the general "C„ interferometer with a two dimensional baseline configuration, 
all the fringes are superposed on a two dimensional detector and they cannot be dispersed. To keep the 
fringe visibilities high, the fractional bandwidth is limited approximately to D/B where D is the aperture 
diameter and B is the length of the longest baseline. In the "C^ interferometer with a one dimensional 
baseline configuration, dispersed fringes are superposed on a two dimensional detector and the fractional 
bandwidth, Au/i', is the same as for the "C2 interferometer. Therefore 




(42) 



For spectrally narrow light, the spectral band widths do not matter and < Ci >=< C2 >. Then 




(43) 



(44) 
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where 6*2 is the total number of photoelectron count for the "C,j interferometer. Since 



ft/ft. ^, 



the SNR ratio of the "C2 and "C„ interferometers is given by 



SNR1/SNR2 = 



'<Ci> n 
< C2 > 2(n - 1) 



For a short basehne (small B/D), the "C,i interferometer is superior and for a long baseline, the "C2 
has an advantage. From SNRi = SNR2, we obtain {B / D) cross over given by 

i^) _^2(n-l) 

V /cross over — * ■ V 'y 

D Av n 

For a continuum source, the fractional bandwidth Aiy/iy for which the source structure does not change 
is the order of 0.1. For ?i = 5, [B / D)cross over — 16. 

From (43), the ratio of SNR for the "C2 and "C,^ interferometers is given by 



SNR,/SNR', = ^i^^;^y (48) 

as in the case of spectrally narrow light. To be fair, we should compare the same one dimensional 
baseline configuration with different beam combination geometries. In this case, the "C^ interferometer is 
superior to the "C2 interferometer by at most a factor ^/2. 

C. Read noise limit 

In this subsection, we evaluate the numbers of detector pixels. Pi, P2 and Pj. There are 2v j Av fringe 
cycles in each intcrferogram. For the Nyquist sampling theorem, '^^v j Av pixels are needed in the fringe 
direction. For the "C2 and "C^ interferometers, we consider the number of pixels in the cross fringe direction 
(dispersion direction) in each dispersed fringe. The fractional spectral bandwidth Av* jv for which bandwidth 
smearing is small, is given by Av* jv < D/B. Therefore for a spectral bandwidth of Av, we need Av j Av* > 
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/ v) / {D / B) pixels in the dispersion direction. Here we take 2{/S.v / v) / [D / B) pixels as the number of 
pixels in the cross fringe direction. So the minimum number of pixels of the two dimensional detector is 
Pi = = 4i^/Ai^ X 2{l\v / v) / {D / B) = 8(B/D). For the "C„ interferometer, the longest baseline can 
be aligned to one of the two sides of a two dimensional detector and the minimum number of pixels is 
P2 = (AiV I t^vY ~ 16(i?/-D)^, and the maximum allowed fractional bandwidth is D/B. 
In the read noise limit, 



SNR, . <C.>^mit)V¥^n 



and 



SNR', = l^k^SllI^M^EIIE, (51) 

where (44) and (45) hold among < Ci >, < C2 >, and < C2 >. For a two dimensional baseline 
configuration, the SNR ratio of the "C2 and "C„ interferometers is given by 

SNRJSNR, = ^^(^)3/2^ (52) 

For a short baseline the "C„ interferometer is superior and for a long baseline the "C2 interferometer 
is better. From the condition that SNRi — SNR2, 

(-) r ("-^)% " )2i^ (53) 

\jjJcross over — y ^ ^ Ai/ J ' 

For n = 5, ba/jv = 0.1, {B / D)cross over ~ 9.3. For a one dimensional baseline configuration, the SNR 
ratio of the "C2 and "C'j interferometers is given by 



SNRi/SNR'2 ^ (54) 
n — 1 
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Since SNRi < SNR'2 for n > 3, the "C'^ interferometer is always superior to the "C2 interferometer 
in the read noise hmit. 

In summary, the "C„ interferometer with a two dimensional baseline configuration has a limited use 
for a short baseline interferometer or possibly a deployable space interferometer. For a separated spacecraft 
interferometer with long baselines, the "C2 is more suitable and the "C^ interferometer with a one dimen- 
sional baseline configuration is even better. In the background limit, the advantage of "C^ interferometer is 
small, while in the read noise limit, that is significant. This is simply because of the difference in the number 
of detectors or that of detector pixels which affect the total read noise. 

5. Sensitivity in the Presence of Natural Thermal Background 
A. Natural Thermal Background 

We are most interested in the potential of an infrared interferometer in the presence of unavoidable natural 
thermal background. Here we consider the emission from interplanetary dust (IPD) and interstellar dust 
(ISD) as the natural thermal background. 

The brightness of IPD emission depends on ecliptic coordinates." Lockman Hole (A = 137°, /3 = 46°) 
which is close to the ecliptic pole is the representative of the dark portion of the sky and the ecliptic plane 
(A = 122°, j3 = 0°) is that of the bright portion of the sky apart from the galactic plane. IPD emission of the 
latter is brighter than that of the former by a factor 3 to 4. The spectrum obtained by the Diffuse Infrared 
Background Experiment (DIRBE) on board the Cosmic Background Explorer (COBE) of the IPD emission 
we use is given in Rcf. 11. 

The brightness of ISD emission is correlated with the distribution of the hydrogen column density. The 
interstellar neutral hydrogen column density of the Lockman Hole is 4—6x10^'^ atoms m^^ and that of the 
ecliptic plane is 4.0x10^"' atoms m~^. To evaluate the dust thermal emission associated with the neutral 
hydrogen columns, we use the brightness per hydrogen atom at 154 fim obtained by a rocket-borne telescope, 
vI^{l^A)/N{HI) — 3.2 X 10"'^^ W sr^^, the dust temperature of 16.4 K, and the gray-body dust emissivity 
law proportional to }^ 

In the calculations below, we use the natural thermal background for the Lockman Hole and note here 
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that the hniiting flux is higher by about a factor of 2 for the ecHptic plane when the detection is Hmited by 
the natural thermal background. 

B. Separated spacecraft interferometer like TPF 

For TPF, a four element array with 3.5 m diameter apertures is considered. The telescope apertures are 
cooled by radiation to 40 K. Here we assume the emissivity of the apertures to be 5%. Varying the length of 
the longest baseline B, we calculate the SNRs of the "C2 interferometer for n — i, D ~ 3.5 m, throughput 
T = 0.1, cr = 2e^, integration time t ~ 100 seconds, and Ai'/i' = 0.1. The read noise of cr = 2e~ may appear 
small, but such a low level after multiple sampling is set as a goal for NGST.^'^ The coherent integration 
time, t, is somewhat unknown, but t much longer than 100 seconds is unlikely due to disturbances on the 
spacecrafts. 

In Fig. 4, 5 fj detection limits of the "C2 (solid line) and "C,'j (dashed line) interferometers for 3600 
second observations are given for baseline lengths, 100 m and 1 km. Read noise is significant at A < 5 fim. 
In the read noise limit, the detection limits of the "C2 interferometer with a 100 m baseline and the "C^ 
interferometer with a 1 km baseline are almost the same. Around 10 /im, the detection limits are set by the 
natural thermal background and at A > 25 /im, they are set by the telescope thermal background. For the 
thermal background limit, the detection limits are baseline independent, because the thermal background is 
independent of the number of detector pixels. It is noted that there will be no far-infrared capability of TPF 
itself according to the current concept. 

C. Separated spacecraft interferometer like Darwin 

For Darwin, a six element array with 1.5 m diameter apertures is considered. Radiation cooling is expected 
to cool the telescope apertures to 30 K. Here we again assume the emissivity of the apertures to be 5%. We 
calculate the SNRs of the "C2 and "C,'j interferometers ior n = 6, D = 1.5 m, T = 0.1, cr = 2e^, integration 
time t = 100 seconds, and Ai^/jy = 0.1. In Fig. 5, 5 cr detection limits for 3600 second observations are given 
as functions of wavelength for baseline lengths, 100 m, and 1 km. For 1 km baseline, the effect of read noise 
is visible up to 10 /.(m. With the number of apertures, n = 6, the "C,'j interferometer with 1 km baseline is 
more sensitive than the "C2 interferometer with 100 m baseline in the read noise limit. 
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D. Single spacecraft interferometer 

Here we consider a deployable interferometer with five 1.5 m apertures and the two dimensional "C„ beam 
combination geometry. We assume actively cooled telescope apertures of 5 K and emissivity of 5%. The 
system throughput T is assumed to be 0.1. We also assume a = 2e~ , integration time t = 100 seconds, and 
Aiy/jy = D/B. We consider the baseline lengths of 10, 30, and 100 m. In Fig. 6, 5 cr detection limits for 3600 
second observations are given as functions of wavelength. The baseline length of 100 m may be too large for 
the "C„ geometry due to the narrow fractional bandwidth and the effect of read noise. The two dimensional 
"C„ interferometer has an advantage for a compact array as discussed above. 

6. Capability for a deep far-infrared survey and comparison with observatories 

In this section, we investigate the potential of a TPF-like space interferometer for a deep far-infrared survey 
which is important in studying the galaxy formation history. 

We first summarize the current status of the subject. In Fig. 7, source confusion limits at 90 and 170 /im 
of different observatories in space predicted by the model E of Guiderdoni et al.^'^ are compared with those 
estimated from the recent Infrared Space Observatory (ISO) observations^^ '.^^ The model E of Guiderdoni 
et al. is known as one of the very successful models which can explain the Infrared Astronomical Satellite 
(IRAS) 60 fiin counts, ISO 15 /im counts, and cosmic far-infrared background radiation recently detected by 
the Far Infrared Absolute Spectrophotometer (FIRAS) on board COBE. 

At 90 /im, the true confusion limit is one order of magnitude greater than the limit estimated by the 
model at an angular resolution of 10" or for a telescope diameter of 2.2 m. The model predictions of the 
confusion limits of SPICA or FIRST at both 90 and 170 jim should be regarded as lower limits. 

In Fig. 8, 5 (T detection limits of a TPF-like "C2 interferometer (4 x 3.5 m) for telescope temperatures, 
5, 10, 20, and 40 K with 5% emissivity are compared with those of NGST, SPICA, and FIRST (Nakagawa et 
al. 1998). The assumed total observing time is 3600 seconds. As for the interferometer, we assume T = 0.1, 
cr = 2e~, Aiy/iy = 0.1, coherent integration time t = 100 seconds, and B = 100 m. In Fig. 9, 5 cr detection 
limits of a TPF-like "C^ interferometer are calculated for the same conditions. 

Before interpreting Figs. 8 and 9, the limitations of the comparison between the SNRs of an interferom- 
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cter and a single telescope must be noted. For the SNR of the interferometer, we use that for a point source. 
In reality there may be multiple sources and some of them may be resolved, in which case the complex 
visibility 7exp(i0) must be obtained by (4) and be inserted into the SNR formula. The sensitivity of the 
observatory is also given for a point source. However a single telescope can observe multiple point sources 
at once, while the interferometer needs to cover sufficient uv points or baseline configurations to recover 
multiple sources. So necessary observing time is likely to be greater for the interferometer. We proceed with 
the comparison with this limitation in mind. 

In the near infrared, the "C^ interferometer is more sensitive than the "C2 interferometer. However the 
spatial resolution of NGST is sufficient and its sensitivity is superior to both interferometers for the purpose 
of galaxy count. In the mid infrared, the sensitivities of the interferometers are comparable to the sensitivity 
of SPICA though spatial resolutions are different. 

The advantage of the interferometers is obvious in the far infrared. In the present concept of TPF, a 
telescope temperature of 40 K is assumed. Even with 40 K telescopes, the interferometers are more sensitive 
than SPICA or FIRST at A > 100 fim. It is clear from the the result of calculation for 5 K telescopes that 
an actively cooled interferometer in the far infrared regardless of the beam combination geometry is ideal for 
detecting faint galaxies which are otherwise undetected due to source confusion. Bright starburst galaxies 
will be detectable at 170 /im out to z = 4 with the actively cooled interferometer. The survey may not only 
identify far- infrared luminous galaxies at high redshift as IRAS did at low redshift, but may even detect 
a galaxy population which is undetectable by NGST in the near infrared due to the weakness of stellar 
emission extincted by dust. 

7. Conclusions 

In this paper, we have studied the sensitivities of space borne infrared interferometers for the purpose of 
general synthesis imaging. Our primary conclusions follow. 

(1) We have derived the expressions for the SNRs of "C2 and "C„ interferometers in the presence of 
source shot noise, thermal background shot noise, and detector read noise. 

(2) We have investigated the tradeoff between the "C2 and "C„ interferometers for a two dimensional 
baseline configuration and found that the "C2 geometry is suitable for a long baseline interferometer, while 
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the "C„ geometry is fitted to a short baschne dcployable interferometer. For a one dimensional baseline 
configuration, we have found that the "C^ interferometer is superior especially in the case of read noise 
limited observations. 

(3) We have presented the detection limits of separated spacecraft interferometers like TPF and Darwin 
in the presence of the natural thermal background as functions of wavelength. The comparison of a TPF-like 
interferometer with NGST. SPICA and FIRST has revealed that at A > 100 ^m, an interferometer is more 
sensitive than SPICA or FIRST. This is true with only radiation cooling of the apertures. 

We thank David Saint- Jacques and the anonymous referee for comments on the manuscript. TN is 
supported by Grant-in- Aid for Scientific Research of the Japanese Ministry of Education. Culture, Sports, 
and Science (No. 10640239). 

A Variance l^[j2(9)] 

The derivation of the variance F[i2((3')] is analogous to that of but more complicated. We utilize the 

results by PK89 whenever possible. 



VMq)] - < ^2{q? > - < ^2{q) 

= X] X^t"^ Re{yij)Re{yki) > - < Re{yij) >< Re{yki) >] cos{ujijq) cos^iOkiq) 

i<j k<l 

-[< Re{y^.j)Im{yki) > - < Re{yij) >< Im{yki) >]coii{uJijq) sin{ujkiq) 
-[< Imiy^j)Re{yki) > - < Im{yij) >< Re{yki) >]sm{LL>ijq) cos{LL>kiq) 
+ [< Im{y.,j)Im{yki) > - < Im{yij) X Im{yki) >] sin(tj„-9) sin(wMg) (55) 

As in §2, the problem of estimating the variance in the image thus reduces to that of estimating three 
types of covariance term: coY[Re{yij),Re{yki)], C0Y[Re{yij), Im{yki)], coY[Im{yij), Im{yki)]. Using (20), we 
obtain the covariance of the real components of a pair of fringe phasors as 

cov[Re{yij),Re{yki)] = y^(< k{p) > -t-g^) cos(ptj,j) cos(ptJfci) 

p 

= ^ Qo > n + < Qo > ^ jgh cos{pujgh + 4>gh) 
P g<h 
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X cosipuij) cos{pujki) (56) 
Every cosine can be written as the sum of exponential functions: 



E exp(zpcjy) + exp(-ipu;y) exp{ipuJki) + exp{-ipujki) 
cos(ptJy)cos(pwfeO = 2^ — ^ ^ 

p p 

- 



2 2 

lcxp{ip{ LL>ij +ujki)) cxp{ip{ujij - ujki)) 

4 4 
p 

exp{-ip{ujij + ujki)) , cxp{-ip{u}ij - ujki)) 



4 4 
= P^. (57) 

In (57), from nonrcdundancy of baselines, ±ujij±ujki 7^ for i ^ k, j ^ I. Similarly to the two frequency 
case, the three frequency case can be computed. 

Following PK89, we now impose an additional condition to ordinary nonrcdundancy in order to simplify 
the calculations. This condition, hereafter referred to as the nonrcdundancy of triangles, concerns three- 
frequency combinations. Specifically, we assume that 

LUgh ± LOij ± Ukl ^ 0, (58) 

unless (gh), (ij), and (kl) form the sides of a triangle. This condition is easily fulfilled for a space 
interferometer with a small number of apertures. With the nonrcdundancy of triangles for i < j and k < I, 



Igh COs{pUJgh + (kgh) COs{pLJij) COs{pUJkl) = ^ Aij^kl (59) 

p g<h 

where the symbol Aij^fc; unless ij and kl baselines form two sides of a triangle, in which case it equals 
1, and 7exp(i0) is the phasor on the third side of that triangle. 
From Eqs (57) and (59), 



cov[i?e(y,,), i?e(yfcO] = ^^^-^^^-^5,^5,1+ < M > 2£^A,,, (60) 

Similarly, 
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C0Y[Im{yij),Im{yki)\ = ^ kkOjiT < M > — - — A.j^ki- (61) 

Everywhere in this section, the upper sign (or expression) is the correct one when the sides (ij) and (kl) 
of the triangle meet at that vertex for which the label has a value intermediate to those of the two vertices, 
i.e., either when i<j = k<lork<l = i<j. The lower sign (or expression) is the correct one otherwise, 
i.e., when i — k or when j = I. 

Evaluation of cov[i?e(?/y ), Im{yki)] is identical with that for the case with Poisson noise alone as obtained 
in PK89: 



coY\Re{yij)Jm{yki)] =< M > ^^^^ A^-fc; x 



where sgn(x) is +1 if a; > and —1 if a; < 0. 
Furthermore, 



+1 

sgn{i — j) for i = k, sgn{i — k) for j = I 



} , (62) 



cov[Re{yij), Im{yki)] ^ ±cov[Im{yij), Re{yki)]. (63) 
Now we evaluate the variance. 



V[i2{q)] = 2^[cos {quj^j) + sm {qoj^j)] 



2 

i<3 



< M > 



1^ ^ Aij^ki cos[g(wy ± ujki)h cos </> 



2 

i<j k<l 

- ^ sin[q{LUij + ujki)hii sin (l)ii 

i<_j=k<l 

- ^ sin[q{ujij + ujki)hkj sin (j)k j 

k<l=i<j 

- ^(n - Z + j - 1) sin[q(ujii + ujij)]'yji sin (j>ji 

3<l 

- ^ i + ^ ~ 1) sii\[q{---ujii + sinipij. (64) 

l<3 

By relabeling the indices slightly and combining the various sums using simple trigonometric identities, 
we obtain the final expression of [12 (<?)], (34). 
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B Remapping of Pupils 

The anonymous referee pointed out the possibihty of pupil remapping which changes the interpretation of 
the interferometer types. We do not consider this possibihty in the main text, but mention here the view 
based on it. 

The distinction between the "C2, "C„, and "C2 refers actuaUy to the "output pupil" (arrangement of 
the beams immediately before fringe detection) rather than to the configuration of the interferometer as a 
whole. One could use a two dimcnsionally distributed array with a linear nonredundant arrangement of the 
beams before the detector like the one shown in Fig. 3, with a remapping stage during data processing. Any 
two dimensional interferometer can therefore accommodate spectral dispersion in principle. 

However such a scheme would require additional optics, affecting the system throughput and thermal 
environment. It also makes the image reconstruction procedure complex. We therefore would not include 
this scheme in this comparative analysis. 
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Fig. 1. Focal plane interferometer and a detector array. This configuration exists for each baseline of an "C2 
interferometer. The effect of shot noise on photodetection process has been fully formulated for this situation. 
For spectrally broad light, the fringe can be dispersed in the cross fringe direction if a two-dimensional 
detector is used. 

Fig. 2. Two-dimensional ^Ca interferometer. Fringes are superposed at the focal plane of a two dimensional 
detector. Since the fringes cannot be dispersed, the spectral bandwidth of this interferometer is limited. 

Fig. 3. One-dimensional '^Cg interferometer. One dimensional fringes are dispersed in the cross fringe direction 
and superposed at the focal plane of a two dimensional detector. Spectrally broad light can be observed by 
an "C„ interferometer with one dimensional baselines. 

Fig. 4. Detection limits of a TPF-like interferometer. 5 a detection limits are plotted as functions of wave- 
length for the maximum baseline lengths given in the figure for an "C2 interferometer (solid) and a one- 
dimensional "C^ interferometer (dash). Assumptions about other parameters are given in the text. 

Fig. 5. Detection limits of a Darwin-like interferometer. 5 cr detection limits are plotted as functions of 
wavelength for the maximum baseline lengths given in the figure for an "C2 interferometer (solid) and a one 
dimensional "C^ interferometer (dash). Assumptions about other parameters are given in the text. 

Fig. 6. Detection limits of a five element compact array with artificially cooled apertures. 5 a detection limits 
of a two-dimensional "C„ interferometer are plotted as functions of wavelength for the maximum baseline 
lengths given in the figure. Assumptions about other parameters are given in the text. 

Fig. 7. Confusion limits of monolithic observatories. The confusion limits predicted from the model E of 
Guiderdoni et al.^"* (solid line) are compared with those estimated from ISO observations by Matsuhara et 
al.^^ (polygon) and Puget et al.^'^ (triangle). 

Fig. 8. Comparison with monolithic observatories. The sensitivity of an "C2 interferometer composed of four 
3.5 m apertures is plotted for telescope temperatures of 5, 10, 20, and 40 K with emissivity of 5%. The 
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maximum baseline length is assumed to be 100 m. Dashed lines denote the sensitivities of NGST. SPICA, 
and FIRST given for comparison. 

Fig. 9. Comparison with monolithic observatories. The sensitivity of a one dimensional "C,j interferometer 
composed of four 3.5 m apertures is plotted for telescope temperatures of 5, 10, 20, and 40 K with an emissivity 
of 5%. The maximum baseline length is assumed to be 100 m. Dashed lines denote the sensitivities of NGST, 
SPICA, and FIRST given for comparison. 
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This figure "Figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/00 11 309v 1 



This figure "Figure3.gif" is available in "gif" format from: 
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